Design, Synthesis and Cytotoxicity Evaluationof New 1,2-diaryl-4, 5, 6, 7-Tetrahydro-1H-benzo[d] Imidazolesas Tubulin Inhibitors.

A new series of 1,2-diaryl-4,5,6,7-tetrahydro-1H-benzo[d]imidazoles, possessing trimethoxyphenyl pharmacophore, were synthesized to evaluate their biological activities as tubulin inhibitors. Cytotoxic activity of the synthesized compounds 7a-f was assessed against several human cancer cell lines, including MCF-7 (breast cancer cell), HEPG2 (liver hepatocellular cells), A549 (adenocarcinomic human alveolar basal epithelial cells), T47D (Human ductal breast epithelial tumor cell line) and fibroblast. According to our results, HEPG2 seems to be the most sensitive, while MCF7 was the most resistant cell line to the compounds. All the compounds expect 7b, possessed satisfactory activity against HEPG2 with mean IC50 values ranging from 15.60 to 43.81 µM.


Introduction
Based on WHO reports, cancer will be the first cause of death in the entire world in future (1,2). There is an intense effort in cancer research to develop novel compounds that are capable to stop cancer cell's growth. Tubulins are rather new targets for cancer therapy, because of their crucial role in cellular functions such as cell shape maintenance, body structure of cilia flagella and immature blood vessels, transport and location of organelles (such as mitochondria, vesicles, etc.) and segregation of chromosome during mitosis (3-7).
Two globular tubulin heterodimer subunits (α-and ß-tubulin) alternately arrange to each other in a head-to-tail fashion to form a string which is called protofilament. Twisting around a hypothetical axis, usually 13 protofilaments constitute a hallow microtubule cylinder (8-13), which as mentioned, has a key role in regulation process.
There are some antimitoic agents which act by targeting three sections of ß-tubulin. These sections, thus, was named based on the interaction of the different classes of antimitoic agents with their binding sites on tubulin. The binding regions discussed are vinca, taxol and colchicine sites (14-16).
Colchicine (Figure 1, 1) is one of the earliest microtubule-targeting agents which were isolated from the meadow saffron colchicium autumnale. In fact tubulin was first purified based on its high affinity binding with colchicine and was named as a colchicine-binding proteins reported that in their structures ring A and B are connected to each other in cis geometry via a tetrahydrobenzimidazole ring ( Figure  1, compounds 7a-f). Docking simulation was performed to position these compounds into colchicine binding site to determine the probable binding model, which suggested the possible inhibition mechanism. Cytotoxic activity of the synthesized compound was also studied against MCF-7, HEPG2, A549, T47D and normal fibroblast cells.

Chemicals
All chemicals, reagents and solvents used in this study were purchased from Merck AG and Aldrich Chemical Company and were used without further purifications. Melting points were determined with a Thomas-Hoover capillary apparatus. Infrared spectra were acquired using a Perkin Elmer Model 1420 spectrometer. A Brucker FT-500 MHz instrument (Brucker Biosciences, Germany) was used to acquire 1HNMR spectra with TMS as internal standard. Chloroform-D and DMSO-D 6 were used as solvents. Coupling constants (J) values are calculated in hertz (Hz) and spin multiplicities are given as s (singlet), d (double), t (triplet) and (17)(18)(19)(20)(21). Because of its high adverse effects, the clinical development of colchicine for cancer treatment has not been successful (22). The therapeutic value of colchicine is well established in the treatment of acute gouty arthritis, familial Mediterranean fever (FMF), Behcet's disease and recurring pericarditis with effusion (21,23,24).
In 1989, Pettit and his co-workers isolated combretastatin (CA-4) from the African willow tree, combretumcaffrum (25). CA-4 ( Figure  1, 2) has a great antitumor effect by binding to colchicine binding site on ß-tubulin. It is capable of rapid and selective disruption of the abnormal tumor vasculature and also exhibits potent cytotoxicity against a board spectrum of cancer lines, including those that are multidrug resistance (MDR) (26)(27)(28). Due to the poor aqueous solubility and bioavailability of CA-4 (29),a wide range of structural analogs of CA-4 have been synthesized and structure activity relationship (SAR) has been determined for these derivatives. According to the existing SAR for CA-4, the 3,4,5-trimethoxyphenyl ring (A ring),cis-configuration of the two phenyl rings, two-atom distance between them and a phenyl ring (B ring)with H-donor groups on positions 3 or 4 are essential for optimum activity (30).
In our study, a series of new analogs are

and Colchicin
The mass spectral measurements were performed on a 6410Agilent triple quadrupole mass spectrometer (LCMS) with an electrospray ionization (ESI) interface.

Anticancer activity
Fibroblast, MCF7 (human breast cancer), HEPG2 (human liver carcinoma), A549 (human lung cancer) and T47D (human breast cancer) cell lines were purchased from Pasteur Institute, Tehran, Iran. The cells were grown in RPMI1640 medium at 37 °C under 5% CO 2 supplemented with 10% heat inactivated fetal bovine serum (FBS), 100 U/mL penicillin and 100 µg/mL streptomycin. Cell viability was evaluated using MTT method which is based on the reduction of 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) dye to purple formazan crystals by mitochondrial succinate dehydrogenase enzyme in living cells. The cells were seeded into 96-well plates at a concentration of 6000 cells/well and allowed to incubate for 24 h. The medium was then discarded and different concentrations of the test compounds dissolved in complete medium, were added to each well. After further incubation for 24 h at 37 °C, the medium was aspirated and 100 µL MTT (2 mg/mL) was added to the wells and incubated for 3 h at 37 °C. The produced formazan crystals were dissolved in 100 µL of DMSO. Plates were incubated for 20 min at 37 °C and the optical densities were read at 570 nm with a reference wavelength of 630 nm as background using a spectrophotometer plate reader (Infinite® M200, TECAN). Colchicine was used as positive control and DMSO as the solvent of the test compounds and its final concentration was less than 0.2%. IC 50 was calculated by calibration curve using Prism software. All the tests were performed in triplicates.

Molecular modeling (docking) studies
All the synthesized compounds were subjected to molecular modeling study and their interactions with the colchicine binding site of β-tubulin were simulated. The tubulin structure was downloaded from the PDB data bank (http://www.rcsb.org/ -PDB code: 1SA0) (31, 32). All the compounds were built using Chem Draw software and minimized subsequently. The protein structure was prepared for docking using AUTODOCK Tool. Docking was performed by AutoDock 4.2 program using the implemented empirical free energy function and the Lamarckian genetic algorithm (LGA)(32). Crystallized ligand was removed from crystal protein (1SA0). Polar hydrogens were added and non-polar hydrogens were merged and finally Kollman united atom charge and atom type parameters were added to 1SA0. Grid map dimensions (126×126×126) were set surrounding active site. Lamarckian genetic search algorithm was employed and docking run was set to 50. The results showed good superimposition of the trimethoxyphenyl moiety of the synthesized compound 7a with the corresponding group of colchicine (Figure 3).
Cytotoxic activity of the synthesized compounds 7a-f was assessed against several human cancer cell lines, including MCF-7 (breast cancer cell), HEPG2 (liver hepatocellular cells), A549 (adenocarcinomic human alveolar basal epithelial cells), T47D (Human ductal breast epithelial tumor cell line) and fibroblast following a previously described procedure (34). The results are shown in Table 1. The comparative cytotoxic activities of the synthesized compounds against Fibroblast, MCF7, HEPG2, A549 and T47D indicate difference in responsiveness/sensitivity of different cancerous cells. In general HEPG2 seems to be the most sensitive, while MCF7 was the most resistant cell line to the compounds. All the compounds expect 7b, possessed satisfactory activity against HEPG2 with mean IC 50 values ranging from 15.60 to 43.81 µM.
Compound 7e exhibited the most potent activities against HEPG2. Selectivity of the cytotoxic activity of the six synthesized compounds was assessed by comparing the cytotoxic activity (IC 50 ) of each compound against each cancerous cell with that of the normal human Fibroblast cells (Table 1). Results were expressed as selectivity index (SI). Compounds 7c, 7d and 7e appear to be the most selective against HEPG2. The highest selectivity was, however, observed for compound 7c against HEPG2 cell line (SI~ 4.1).
Molecular modeling studies comparing colchicine and the synthesized compound  7a were carried out. The structures were superimposed by a least square fit of the 3,4,5-trimethoxyphenyl moieties. The results are shown in Figure 3. As anticipated, the crucial orientation and spacing of the two aromatic rings were essentially retained in the new compounds. The bridge moieties cause the A ring and B ring of these compounds to lie in a similar positions to that of the colchicine site. (13)